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ABSTRACT 

We present the first simultaneous infrared (IR) and visible spectro-polarimetric observations of a solar internet- 
work region. The Fe I lines at A6301.5 A, A6302.5 A, A15648 A, and A15652 A were observed, with a lag of 
only 1 min, using highly sensitive spectro-polarimeters operated in two different telescopes (VTT and THEMIS 
at the Observatorio del Teide). Some 30% of the observed region shows IR and visible Stokes V signals above 
noise. These polarization signals indicate the presence of kG magnetic field strengths (traced by the visible lines) 
co-existing with sub-kG fields (traced by the infrared lines). In addition, one quarter of the pixels with signal 
have visible and IR Stokes V profiles with opposite polarity. We estimate the probability density function of 
finding each longitudinal magnetic field strength in the region. It has a tail of kG field strengths that accounts for 
most of the (unsigned) magnetic flux of the region. 

Subject headings: Sun: magnetic fields - Sun: photosphere 
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1. INTRODUCTION 

Most of the solar surface appears non-magnetic when it is ob- 
served in routine synoptic magnetograms (e.g., those obtained 
at Kitt Peak' or with the MDI insti'ument on boai'd of SOHO^). 
However, magnetic fields are detected almost everywhere, also 
in internetwork (IN) regions, when the polarimetric sensitivity 
and the angular resolution exceed a threshold. These magnetic 
fields of the quiet Sun, elusive but now accessible to many ex- 
isting spectro-polarimeters, have received an increasing interest 
during the last years. They produce only weak signals, but they 
cover much of the solar surface and, therefore, they may be car- 
rying most of the unsigned magnetic flux and energy existing 
on the solar surface at any given time (e.g., Sanchez Almeida 
2003, and references therein). This makes the IN potentially 
important to understand the global magnetic properties of the 
Sun. 

The observational study of the IN fields began in the seven- 
ties (Livingston & Harvey 1975; Smithson 1975), but it is still 
in an initial phase due to the complexity of the magnetic topol- 
ogy. Opposite polarities seem to coexist in one and the same 
resolution element of typically l"(see Sanchez Almeida et al. 
1996; Sanchez Almeida & Lites 2000; Sigwarth et al. 1999; 
Lites 2002; Khomenko et al. 2003). In addition, large system- 
atic differences of intrinsic field strength are deduced depend- 
ing on the diagnostic technique. Techniques based on visible 
lines favor kG field strengths (Grossmann-Doerth et al. 1996; 
Sigwarth et al. 1999; Sanchez Almeida & Lites 2000; Socas- 
Navarro & Sanchez Almeida 2002; Dominguez Cerdena et al. 
2003a,b), whereas inferences based on infrared (IR) lines indi- 
cate sub-kG field strengths (Lin & Rimmele 1999; Khomenko 
et al. 2003). (We add for completeness that estimates based 
on the Hanle effect indicate the presence of even weaker fields, 
e.g., Faurobert-Schofl et al. 1995; Bianda et al. 1999.) Sanchez 
Almeida & Lites (2000) conjectured that this systematic dif- 
ference actually reflects the existence of a continuous distri- 
bution of field strengths ranging from sub-kG to kG (see also 

' http://www.noao.edu/kpno/ 

^ http://sohowww.nascom.nasa.gov/ 



Socas-Navarro & Sanchez Almeida 2002). Different observa- 
tions are biased towards a particular part of such a distribution, 
producing the observed systematic difference of field strengths. 
Even more, Sanchez Almeida & Lites (2000) put forward a 
specific physical mechanism that biases the IR observations to- 
wards sub-kG fields. It is caused by the vertical gradient of field 
strength existing in any magnetic structure (provided they sat- 
isfy some sort of mechanical balance with the mean solar pho- 
tosphere). The Zeeman splitting of the IR lines is so large that 
the absorption in various layers, with different field strengths, 
occurs at substantially different wavelengths. The absorption 
is thus spread over a wide range of wavelengths, which dilutes 
the polarization signals. On the other hand, sub-kG fields do 
not disperse the IR absorption. In this case the IR lines are in 
the so-called weak field regime so, independently of the field 
strength, the absorption is always produced at the same wave- 
lengths, and a significant polarization signal builds up. Socas- 
Navarro & Sanchez Almeida (2003) show that the existence of a 
horizontal gradient can also explain the observed discrepancy. 
In this particular case a distribution of field strengths is chosen 
such that the visible lines reflect kG fields because most of the 
magnetic flux is in the form of kG fields. In the IR lines, how- 
ever, the kG signal is much reduced because of the wavelength 
smearing effect described above. 

So far the seemingly contradictory visible and IR observa- 
tions were taken independently (by different authors at differ- 
ent times). If the conjecture on the coexistence of weak and 
strong fields were correct, a simultaneous observation of the 
visible and IR lines would have to give co-spatial signals show- 
ing both, signatures of kG (visible) and of sub-kG (IR) fields. 
The present study was meant to test this specific prediction. 

The simultaneous observations and few details of the data re- 
duction are described in § 2. Then we explain the data analysis, 
the classification of the polarimetric profiles, and the determi- 
nation of magnetic field strengths (§3). We present the main 
results in § 4, and the conclusions in § 5. 
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spectral ranges for each pixel of our FOV. Here we only use in- 
tensity and circular polarization. The noise in the polarization 
measurements was estimated as 3 - 4 x lO""* /c in the IR data, 
and 7 x 10^/c in the visible data (/c is the continuum intensity). 
They correspond to the rms fluctuations of the polarization in 
continuum wavelengths. 

3. DATA ANALYSIS 

We select for analysis those IR and visible Stokes V pro- 
files whose peak polarization is larger than three times the noise 
level. This thresholding separates clear signals from very noisy 
profiles or pure noise. The pixels with both IR and visible sig- 
nals above their respective threshold correspond to 30% of the 
FOV. However, an additional 40% of the pixels present polar- 
ization signals only in either the visible or the IR. This study 
deals with those pixels with signals in both spectral ranges 
(some 1840). 

All selected profiles were classified using a Principal Com- 
ponent Analysis (PCA) algorithm, similar to the one used by 
Sanchez Almeida & Lites (2000, § 3.2). We classify, simul- 
taneously, the IR and visible Stokes V profiles. All vectors are 
scaled to the maximum of the blue lobe of the Fe I A6302.5 A 
Stokes V profile. The PCA classification results in 25 different 
classes with more than 0.5% of the classified pixels, and 15 
classes containing rare profiles (between 0.5% and 0.2%). Fig- 
ure 2 shows some of the significant classes. We note, firstly, 
that the IR profiles in class are similar to the most abundant 
class obtained by Khomenko et al. (2003, Fig. 1). Secondly, 
class 2 corresponds to network structures in our FOV, and it has 
characteristic broad IR profiles (also observed in plage regions 
by Stenflo et al. 1987). Finally, the most conspicuous fact is 
that some classes (e.g., # 4 and # 8 in Fig. 2) show opposite po- 
larity in the visible and in the IR lines (the lobes have opposite 
sign). 

In order to estimate magnetic field strengths, we performed 
Milne-Eddington (ME) inversions of the profiles of all individ- 
ual pixels, as well as the mean profiles among those in each 
class. We use MILK, a code which is part of the Commu- 
nity Inversion Codes'' developed by HAO (see Socas-Navarro 
2001). Each inversion provides 11 parameters, including the 
longitudinal component of the magnetic field strength, and the 
fill factor (fraction of resolution element covered by magnetic 
fields). (For a description of the procedure and approximations, 
see, e.g., Skumanich & Lites 1987.) The inversion is separately 
done for the two pairs of lines, i.e., we carry out two indepen- 
dent inversions per pixel, one for the visible pair and another for 
the IR pair. In each inversion we use the same atmospheric and 
magnetic field parameters for the two lines, whereas the ratio of 
absorption coefficients is fixed according to the expected val- 
ues. Note that our ME inversion does not fit line asymmetries, 
although they are clearly present among the observed profiles. 

4. RESULTS 

The PCA classification indicates that 25% of the pixels with 
coincident visible and IR signals give visible and IR profiles 
with opposite polarity. It follows that the visible and IR fines 
are tracing different magnetic structures in the same resolution 
element. Some classes have IR profiles with three lobes while 

^ Actually, we allow for a shift, a rotation of the FOVs, and a scale factor in the direction along the slit. The rotation is only 1°, whereas the scale factor is consistent 
with the ratio of effective focal lengths in the focal plane of the two polarimeters. 
http://www.hao.ucar.edu/pubUc/research/cic/ 



2. OBSERVATIONS 

We took advantage of the Spanish Observatorio del Teide 
(Tenerife, Spain) having two leading telescopes with sensi- 
tive spectro-polarimeters. The IR observations were carried 
out with the Tenerife Infrared Polarimeter (TIP, Martinez Pil- 
let et al. 1999) at the German Vacuum Tower Telescope (VTT). 
We obtained spectra of all four Stokes parameters for the IR 
lines Fe I A15648 A (effective Lande factor, geff = 3) and 
Fe I A15652 A (geff = 1.53). The observations in the visi- 
ble range were gathered with the spectro-polarimetric mode 
of the THEMIS telescope (Lopez Ariste et al. 2000). We 
chose 3 different spectral ranges for the visible observations. 
Only one containing the lines Fe I A6301.5 A (g^ff = 1.67) and 
Fe I A6302.5 A (geff = 2.5) is analyzed here. The sht width was 
set to 0."5 in both telescopes, and we scanned the solar surface 
from east to west with a step size of 0"5. Each scan consisted 
of 60 positions, with an integration time of 30 sec per position. 

The dataset for this work was obtained on August 10, 2002, 
in a quiet Sun region close to disk center (/j. = 0.94). We se- 
lected IN regions trying to avoid bright features in the Ha and 
Ca n K sUt-jaw filtergrams suppUed by the VTT. Having a few 
network patches in our fairly large Field-of-View (FOV) is un- 
avoidable, though. In order to point both telescopes to the same 
region, we benefited from a video link that provided VTT slit- 
jaw images at THEMIS. First, we pointed both telescopes to an 
easily identifiable structure, e.g., a sunspot. Then we moved the 
VTT pointing to a quiet IN region and displaced THEMIS by 
the same distance. For the data set analyzed here, this method 
of differential pointing gave an offset of the co-alignment of the 
two telescopes of 1". When the offset is corrected it leads to 
a time lag between the IR and visible data of only 1 minute. 
We found this offset by shifting the continuum intensity images 
obtained from the visible and the IR scans (Fig. 1). The op- 
timum shift is chosen as that yielding the minimum difference 
between the two images.-' Figure 1 shows the resulting intensi- 
ties in the IR (left) and the visible (right). Despite the difference 
in angular resolution, it is clear that both images correspond to 
the same region. The white boxes point out a set locations on 
the Sun, and one can easily identify the same structures in the 
boxes of the two spectral ranges. The spatial resolution for the 
IR observations is 1"-1 ."2, with a granulation contrast of 1 .7%, 
and it becomes 1 ."5-1 "7 in the visible, with a contrast of 1 .9% . 
The resolution was estimated as the cutoff of the power spectra 
of the continuum images shown in Fig. 1 . The IR resolution 
is better due to the use of an active tip-tilt mirror that removes 
image motion during scaiming (Ballesteros et al. 1996). The 
FOV connmon to visible and IR observations is of the other of 
30" X 35" (Fig. 1). 

The data reduction for the IR and visible ranges was very 
similar: flat field correction, removal of fringes, demodula- 
tion to get the Stokes parameters from the measurements, and 
correction of seeing-induced crosstalk by combining the two 
beams of each polarimeter. The instrumental polarization (IP) 
of the VTT was cleared away using a model Mueller matrix for 
the telescope, which was calibrated with large linear polarizers 
inserted at the telescope entrance (CoUados 2003). THEMIS 
data do not need this step since it is an IP free telescope. The 
reduction process renders Stokes /, Q, U and V profiles in both 
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the visible profiles have only two (not included in Fig. 2). Ac- 
tually, the observed Stokes V profiles present a wide variety 
of asymmetries, in agreement with previous observations (e.g., 
Sanchez Almeida & Lites 2000; Khomenko et al. 2003). 

The ME inversions show that most of the classes have visi- 
ble Stokes V profiles characteristic of kG fields, while, simul- 
taneously, the IR profiles indicate the presence of sub-kG field 
strengths. This also happens when we analyze all individual 
profiles. The result cannot be significantly biased by noise, 
since the mean profile of a class results from averaging hun- 
dreds of profiles and, therefore, bears negligible noise (see Fig. 
2). On the other hand, the systematic difference between the 
IR and the visible is insensitive to the specific technique used 
to diagnose the field strength. The ratio between the Stokes V 
extrema of the two visible lines can be used to diagnose the 
regime of field strength in the atmosphere. A ratio close to 1 
indicates kG magnetic fields, whereas a value of some 0.6 is 
expected for sub-kG fields (see Sect. 5 in Socas-Navarro & 
Sanchez Almeida 2002). Most of the classes have ratios be- 
tween 0.9 and 1, which evidences IN magnetic structures in the 
regime of kG fields. After the inversion of profiles in the vari- 
ous classes, and also from individual profiles, we found a mean 
magnetic field strength of 1 100 G from the visible lines. They 
fill only 1% of the resolution elements. The mean field strength 
in the IR is 300 G, and covers some 2% of the part of FOV 
studied here. 

The ME inversions also show a few classes with weak mag- 
netic fields both in visible and the IR. There are also classes 
with strong fields in both spectral ranges corresponding to net- 
work points in our FOV (class 2 in Fig. 2). 

Longitudinal magnetograms measure flux density. One can 
infer it from ME inversions as the product of the longitudinal 
component of the field strength times the fill factor. When aver- 
aged over the 30% of the FOV that we study, the mean unsigned 
flux density turns out to be 1 1 G in the visible, and 6 G in the IR. 
Due to the complexity of the IN magnetic fields, the observed 
unsigned flux is expected to increase with increasing angular 
resolution (see, e.g., Sanchez Almeida et al. 2003; Domihguez 
Cerdena et al. 2003b). The fact that the IR, with higher reso- 
lution, has lower unsigned flux, seems to indicate that this is a 
true solar effect. Specifically, most of the unsigned magnetic 
flux of the IN seems to be in the form of kG fields. On the 
other hand, the signed flux of the region in the visible is -1 .5 G, 
while it is less than 0.2 G in the IR. A similar net flux has been 
observed before in visible lines (e.g. Lites 2002; Domihguez 
Cerdeiia et al. 2003b). We find it to be co-spatial with very low 
signed flux in IR fines. 

Using magnetic fields from the ME inversion, we have es- 
timated the Probability Density Function (PDF) of finding a 
given longitudinal magnetic field strength in our IN quiet Sun 
region. We assume that the visible and the IR polarization 
signals trace different structures (a reasonable hypothesis, in 
view of the findings described above). The range of observed 
field strengths is divided in bins of 50 G. Then, treating visi- 
ble and IR observations as independent measurements, we as- 
signed each measured longitudinal magnetic field to one of the 
bins. The probability of finding magnetic field strengths in 
a given bin (i.e., the PDF) will be proportional to the sum of 
fill factors of all those measurements in the bin. The result is 
shown in Fig. 3, the sohd line. This PDF has been normal- 
ized so that its integral is 0.3, i.e., the fraction of FOV under 
analysis. In other words, the normahzation assumes that the re- 



maining 70% of the surface is field free. Figure 3 includes the 
PDF for the IR data alone (the dashed line), and for the visible 
data alone (the dotted line). We also plot the PDF suggested 
by Socas-Navarro & Sanchez Almeida (2003, the dash-dotted 
line). We find that 75% of the total flux is concentrated in fields 
stronger than 500 G, which, however, occupy only 0.4% of the 
full FOV. 

A final comment is in order. The visible and IR maps that we 
employ have a time lag of one minute and were coaligned with 
a finite accuracy (however better than the angular resolution). 
The lack of strict co-ahgnment and simultaneity may produce 
some differences between the visible and the IR spectra, how- 
ever, they would be of random nature and, consequently, they 
can hardly explain the observed systematic differences. 

5. CONCLUSIONS 

The (almost) simultaneous visible and IR Stokes V spectra 
described above have several conspicuous properties. They 
reveal the coexistence of kG and sub-kG magnetic fields 
in most of the pixels with polarization signals above the 
noise (30% of the FOV). The IR lines indicate sub-kG fields 
whereas visible fines trace kG. This demonstrates the conjec- 
ture by Sanchez Almeida & Lites (2000) and Socas-Navarro & 
Sanchez Almeida (2002, 2003). 

At least 25% of pixels with coincident visible and IR signals 
have one polarity in the visible and the opposite in the IR. This 
second result offers a direct proof of the existence of two polar- 
ities in 1 ."5 resolution elements. They have been inferred from 
line asymmetries (Sanchez Almeida & Lites 2000), as well as 
from the skewness of the distribution of line-of-sight magnetic 
field inclinations in measurements away of the solar disk center 
(Lites 2002). 

The Probability Density Function for the longitudinal mag- 
netic fields, that we infer from ME inversions, has a long tail 
of kG fields. Some 75% of the observed (unsigned) magnetic 
flux is in field strength larger than 500 G. This excess of kG 
flux is in agreement with the PDF suggested by Socas-Navarro 
& Sanchez Almeida (2003) to reconcile discrepancies between 
field strengths obtained with different methods (see § 1). We 
have also found the extended tails in the IR line Fe I A15648 A 
predicted by such PDF (see classes and 5 in Fig. 2). They 
carry the information on the strong fields, but are easily misin- 
terpreted due to the large lobes formed by weak fields. 

Our observations prove that the use of only visible lines or 
only IR fines to study the IN fields biases the results. Diagnos- 
tics based on simultaneous visible and IR observations offer a 
viable and refiable alternative. 
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Fig. 1 . — Continuum intensity images from the IR (left) and the visible spectro-polarimetric scans (right). The squares point out structures which can easily be 
identified in both images demonstrating the good co-alignment. We scanned along the horizontal direction (solar East- West); the vertical direction is the direction 
along the spectrograph slit (North-South). 
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Fig. 2.— Representative Stokes V profiles from the PCA classification. Mean profiles of Fe I A6301.5 A (left), Fe I A6302.5 A, Fe I A15648 A, and 
Fe I A 15652 A (right) among those in classes 0, 2, 4, 5 and 8 (see the labels on top of the figures, where the percentage of profiles included in each class is 
also given). Wavelengths are in A relative to the line center (indicated in each figure as a large plus sign). The small circles represent actual observations whereas 
the solid lines correspond to ME fits to the observed data. 
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Fig. 3. — Probability Density Function (PDF) for observing a longitudinal magnetic field strength Bz. The solid line includes both IR and visible data, whereas 
the distribution of values obtained solely from the IR lines and the visible lines are shown as the dashed and the dotted lines, respectively. The dot-dashed PDF was 
suggested by Socas-Navarro & Sanchez Almeida (2003) and it is represented for reference. 



